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We investigate the relation between the BPS baby Skyrme model and its vector 
meson formulation, where the baby Skyrme term is replaced by a coupling between 
the topological current and the vector meson field u^. The vector model still 
possesses infinitely many symmetries leading to infinitely many conserved currents 
which stand behind its solvability. It turns out that the similarities and differences 
of the two models depend strongly on the specific form of the potential. We find, for 
instance, that compactons (which exist in the BPS baby Skyrme model) disappear 
from the spectrum of solutions of the vector counterpart. Specifically, for the vector 
^-i model with the old baby Skyrme potential we find that it has compacton solutions 

only provided that a delta function source term effectively screening the topological 
charge is inserted at the compacton boundary. For the old baby Skyrme potential 
squared we find that the vector model supports exponentially localized solitons, 
like the BPS baby Skyrme model. These solitons, however, saturate a BPS bound 
which is a nonlinear function of the topological charge and, as a consequence, higher 
solitons are unstable w.r.t. decay into smaller ones, which is at variance with the 
more conventional situation (a linear BPS bound and stable solitons) in the BPS 
baby Skyrme model. 

PACS numbers: 11.30.Pb, 11.27.+d 

I. INTRODUCTION 



In the low energy limit QCD becomes a strongly interacting quantum field theory where 
the primary quark and gluon fields are no longer proper degrees of freedom. Indeed, due 
to confinement they form colorless bound states, i.e., mesons and baryons (and perhaps 
glueballs). Although the complete description of these effective particles from first prin- 
ciples, that is by direct computations in QCD, is a very complicated and still unsolved 
problem, several effective theories have been proposed. One of the best known and success- 
ful approaches is the Skyrme framework in which pions are identified as the proper low 
energy degrees of freedom. Then, baryons appear as collective, topological excitations of 
the mesonic, matrix-valued field U G SU(2) with an identification between the baryon num- 
ber and the topological charge Q G n 3 (S 3 ). Unfortunately, the exact form of the effective 
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action is unknown. Typically, one includes the standard sigma model part (kinetic term), 
the Skyrme part (the fourth derivative term which is required by the Derrick theorem) and 
a potential (chiral symmetry breaking term). A systematic derivation of the model from 
QCD has not been completed yet. Under additional assumptions (1-loop expansion, large 
N c limit and small derivative limit) one can derive Skyrme-like models, but with additional 
terms which seem to be rather problematic. They contain higher powers and higher or- 
ders of time derivatives that lead to problems with the (semi-classical) quantization and 
dynamical properties of solutions. Recently, a new Skyrme-type action has been proposed, 
the so-called BPS Skyrme model 0]. It consists of two parts: the baryon density cur- 
rent squared (the sextic term) and the potential. This model possesses many interesting 
mathematical properties, like infinitely many symmetries with infinitely many conserved 
charges, solutions saturating a BPS bound which are given by exact formulas for any value 
of the topological charge. In addition, it cures some of the phenomenological problems of 
the standard Skyrme model (strongly non-BPS solutions, large binding energies, lack of the 
incompressible fluid property). 

A natural questions one may ask, once a Skyrme type model is specified, is how the proper- 
ties of the solitons (baryons) are changed by interactions with the lightest fields i.e., the U(l) 
Maxwell field and the vector mesons w M . Specifically, the interaction between the Skyrme 
field and the vector mesons is given by the coupling of the vector field with the topological 
current density u^B^ Such a coupling is known to effectively induce a term propor- 

tional to the square of the baryon current. Because of the high nonlinearity of the vector 
Skyrme model we do not know much about its solitonic solutions. Usually one is left with 
numerical solutions [7] , and a detailed understanding of the relation between solitons in the 
Skyrme and the vector Skyrme model is still rather incomplete. 

In this work we want to analytically investigate this issue in (2 + 1) dimensions using the 
ow dimensional version of the BPS Skyrme model [s| i.e., the BPS baby Skyrme model 9|- 



In this lower- dimensional case, the term induced by the vector meson (the topological 



current squared) is, at the same time, the usual Skyrme term which is of fourth power in 
derivatives. In particular, we will show for which potentials both models (the BPS baby and 
its vector version) have qualitatively similar solitonic solutions, and for which potentials the 
properties of solutions are rather distinct. 

An additional, related aim of this paper is to shed further light on possible applications 
of the BPS Skyrme model to baryon physics. The BPS Skyrme model with the standard 
potential (i.e., the pion mass term) may be viewed as a chiral theory where the pions are 
very (infinitely) massive in the sense that the normal linear pion excitations are completely 
suppressed. On the other hand, pions are the lightest effective particles. Hence, the inclusion 
of the vector mesons with finite mass and a standard kinetic term in the BPS Skyrme model 
in some sense reverses the typical mass hierarchy and, therefore, it would be surprising if 
it lead to qualitatively similar solitons as in the original BPS model. Here, we show that 
in 2+1 dimensions the solitons in the baby BPS model with the standard mass term for 
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the "pions" (the so-called old baby Skyrme potential), and, therefore, with infinitely heavy 
"baby" pions in the sense just explained, are in fact completely different from the solitons 
in its vector version (with finite mass "baby" vector mesons). These findings should be 
contrasted with an analogous investigation for the full baby Skyrme model 12] with the 
"old" baby Skyrme potential, where the authors find that the solitons in the original baby 
Skyrme model and its vector version behave rather similarly 

Finally, let us mention that one way to pass from the BPS baby Skyrme model to the vector 
model is by replacing the local integration kernel in the Skyrme term with a nonlocal integral 
kernel of the Yukawa or Coulomb type. That is to say, instead of the Skyrme term 

L Sk = Bl(x) = J d 3 yB^x)g^6 3 (x - y)B v (y) 

we have the nonlocal term 

Kec = J d 3 yB^x)K^(x - y)B v (y) 

where K^ v (x — y) is the integral kernel of the Yukawa or Coulomb type which results from 
integrating out the (massive or massless) vector mesons. The same nonlocal interaction 
induced by a Yukawa or Coulomb integral kernel has been investigated for the nonlinear 



Schroedinger equation in 13] 



II. THE u- VECTOR MODEL 

It is well-known known that the coupling with a vector field offers an alternative way 



to stabilize baby skyrmions 12j. Then, instead of the standard (2+1) dimensional Skyrme 
term one introduces a coupling between the topological field and vector mesons by means 
of the topological current B^. Concretely, the Lagrange density reads 

7Y1 2 -» -» I 1 

£ = "y (<V>) 2 - /iV(0) - -{d^ v - d u u,f + -W 1 ^ + AVJ", (1) 
where the topological current is 

B^ = -^-e^$-(dJxd^), (2) 

and <p = ((f) 1 , <p 2 , (p 3 ) is a three component unit vector field. It introduces a nontrivial 
topology into the model since, for static configurations, it may be viewed as a map from 
R 2 U {oo} = S 2 (obtained by a compactification of the base space via adding a point at 
infinity, because <fi — > <p vac when x — > oo) into the target space S 2 . Therefore, we will call 
this field topological, in contrast to the mesonic fields w M . The potential V is a positive 
definite function which depends only on the topological field, and whose particular form 
influences rather strongly the properties of baby skyrmions [l5|. Here, we restrict 

ourselves to potentials which leave some of the global 0(3) symmetry unbroken, i.e., which 
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depend only on the third component of the topological field V = V(<p 3 ). More concretely, 
we will analyze a family of one-vacuum potentials which provides a natural generalization 
of the old baby potential 

where a > 1. 

The extreme (or BPS) version of the model is defined by the same assumption as in the case 
of the usual baby Skyrme model, that is m = 0, i.e., when the standard kinetic term for the 
topological field is absent 

C = -fSVfa) - ~(0„w„ - d u u,) 2 + l -M 2 ul + X'u^ (4) 

Observe that we leave the mesonic part of the model untouched, so it is by no means obvious 
whether all the nice features of the BPS baby Skyrme model will survive also in its vector 
meson formulation. 

For reasons of convenience we prefer to deal with the CP 1 version of the model, that is, we 
apply the stereographic projection 

0= , n ,0 (u + u, -i(u - u), \u\ 2 - l) . (5) 
1 + \ur K ' 



Then, 



C = -^V{\u\ 2 ) - - 4 (d,co u - d u co,r + l -M 2 ul + zAa;^ (i ^ 2)2 (6) 
with potential 

The pertinent field equations are 

dfa + M 2 u u + iXe^ -^L- = 0, (8) 

p (i + \u\ z y 

a-l 

u = 0, (9) 



( 


u\ 


2 \ 


\ 1 + 


\u\ 





and their static versions are 

di u - M u = iX : v , 10 

(l + \u\ z y 

( \u\ 2 \ a_1 _ 

iXCVrtoVchU — V r wV 'au) + a/j 2 — - u = (11) 

\i + \uY j 

where we have used that Bi = 0, Ui = and oj = ojq. An obvious ansatz which is compatible 
with these equations is the axially symmetric one 

u = /(r)e^, u = w(r), (12) 
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where n is an integer number which coincides with the winding number of the topological 
field. Further, the relevant topologically nontrivial boundary conditions are 

f(r = 0) = oo, f(r = R) = (13) 

where, depending on the potential, R may be finite (for compactons) or infinite (for exponen- 
tially or power- like localized baby skyrmions). Additionally, the regularity of the equation 
(finiteness of the energy) for the meson field requires 

co' r (r = 0) = 0, co(r = R) = 0. (14) 

Finally, one usually imposes the first derivative continuity (or regularity for non-compacton 
case) condition 

f( r = R) = 0, u'(r = R) = 0. (15) 

As we will see below, however, in some cases this condition cannot be satisfied. This results 
in an important qualitative difference between solutions of the BPS Skyrme model with the 
old baby potential and its u meson version. 



A. Massless case 

Let us start with the simplest case of massless u mesons. The limit of massless vector 
mesons is not directly related to the induction of the Skyrme term, but it is simpler and 
some features of the solutions will be useful later on. The massless vector meson field may 
also be interpreted as an abelian £7(1) field which here is coupled to the BPS baby model 
by the topological current only. Such a coupling, alternative to the minimal one, does not 
require an unbroken U(l) symmetry for the topological field <fi. Hence, in principle, one 
may consider potentials which completely break the global SO (3) symmetry. Notice that 
a similar non-minimal coupling appears as part of the (3+1) dimensional Skyrme model 



coupled properly to the Maxwell field [16j . Let us also emphasize that the analogy with the 
electromagnetic interaction allows us to conclude immediately that static solutions have to 
obey one of the two following conditions. Either the energy of the static solution is infinite, 
because the electrostatic energy in two space dimensions is infinite (has a logarithmic IR 
divergence). Or the total electric charge and, therefore, the topological charge (which is 
equal to the electric charge in our case) is, in fact, zero, despite our topologically nontrivial 
ansatz. This second possibility implies that the corresponding static solution is really a 
solution in the presence of an inhomogeneous source term forced upon us by the consistency 
of the solution, such that the source effectively screens the topological charge. We shall find 
both possibilities in the sequel. The static field equations are 

1 / \ 2nA ff r 

-d r (ru r ) = -— (16) 
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/ 



r r nX \l + f 2 



a-l 



0. 



Now we introduce the new variable x = r 2 /2 and get 

fSL 



d x (2xu' x ) =2nX- 
where C is an integration constant and 



d x 



2xui' + n\ 



1 + f 



+ C 



f 



/ 2 



a-l 



0. 



(17) 



(18) 



(19) 



The last equation allows to express u' x in terms of the profile function /. Hence, after 
inserting it into (jTSJ) we get a simple algebraic formula for / 



;i + / 2 



a-l 



n 2 X 2 1 X n 



1 + f* 



2a[i 2 x ax 



(20) 



where we have used the boundary condition f(x = 0) = oo to fix the integration constant 
C = 0. Then, the solution should be inserted back into (TT91) to obtain uj. 



1. Compacton - a = 1 
As a first example let us consider the old baby potential i.e., a — 1. Then one finds 



where 



<f - 1 x<X 
fir) { (21) 

x > X 

.d.D: ^ (22) 

x > X 
n 2 \ 2 

X - ^ (23) 

One can see that this solution is only of C-class (i.e. a continuous function without continuous 
derivatives), and, therefore, may introduce a source term at the compacton boundary. 

Proposition 1: Any compacton solution of the massless e.o.m. (lisp . (1191) is a solution 
with a Dirac delta source located at the boundary of the compacton which screens the total 
baryon charge inside the compacton. 

Proof: We start by assuming a compacton solution with a finite compacton radius R 

f(x) x<X 

fir) { (24) 
x > X 
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Cj(x) x < X 

u(x) = { (25) 
x > X 

where f,u solve f|T8|) . f|T9|) inside the ball with the radius X. Specifically, for the massless 
vector model with the old baby potential they are given by fT2~T]) . fl23|) and X = X . Now, 
consider the equation for the meson field inside the compact ball (we use that C = 0) 

d x (W„) = -nA<9, (j^/i) • (26) 

Here, the right hand side, that is the topological charge density, may be treated as a source 
for the mesonic field 

p(x) = —nd x =- ) . 

U + /V 

We easily find the total charge inside the compacton ball 



Qinside = 2tt [ dxp(x) = —2im [ dxd x ( — \ = —Inn 

Jo Hy 1 Jo \ 1 + pj 1 + p 



-27m. (27) 



However, when we look at equation fT26|) at the boundary (or equivalently in the whole space) 
then an additional current appears. Namely, we define a function H(x) 

H(x) = 2xu' x + n ^Y~]~p ( 28 ) 

whose derivative is equal to the static equation for the mesonic field <HM . Due to the 
equations of motion, it holds that H(x) = inside the ball. On the other hand, for the 
vacuum values / = and uj = which the fields take outside the compacton we get 
H{x) = A. Therefore, the function H computed for compacton solutions is proportional to 
a Heaviside theta function 

H(x)\ onshell = n\6(X-x). (29) 

Obviously, its derivative is proportional to the Dirac delta, and a source term at the com- 
pacton boundary emerges. Thus, finally, the compacton solution is a solution to the following 
equation 

d x (2xu' x ) = -n\d x + n\5(X - x). (30) 

The topological charge stored in the compact ball is screened by the source term added at 
the boundary of the compacton. Hence, the total charge of the solution is zero which results 
in the finiteness of the total energy. 
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2. Non-compact solutions - a > 1 



Now, from fl20|) one sees that the profile function tends to at spatial infinity in a power- 
like fashion. Indeed, solving ( 1201) for large x we get as the leading term 



On the other hand, from (|T8l) 



nX 

~2x + "' 



X n \ 2(Q-1) 



x 



+ .. 



nA 

In x + ... 

2 



(31) 



(32) 



and it is precisely this logarithmic behaviour which will give rise to the IR divergent energy. 
To see it, let us consider the static energy integral 



E = 2n I rdr 
Jo 

rod 

2tt / dx 
Jo 



2 T/ /yx 1 2 2An uff' r 



(33) 
(34) 



As always for a " Coulomb" energy, the second and third term combine into a total derivative 
plus a term which would normally give the Coulomb energy. The difference is that here the 
total derivative is IR divergent. Indeed, using the (non-dynamical) field equation (JIHjl (with 
C = 0) for xu' x we get 



E = 2vr / dx 
Jo 



UJ 



nX d 
uj- 



2 dx \l + P 2 dxl + / 2 



(35) 



where the third term is the normal Coulomb term, whereas the second, total derivative term 
is IR divergent in our case due to the u ~ lnx behaviour of oj for large x. 
As an exactly solvable example we may consider the case a = 2. Then 



ft \ X° 

/(«) = Vt' 

c(x) = -^ln(x + X ). 



(36) 
(37) 



Undoubtedly, the massless vector meson version of the BPS baby Skyrme model differs 
significantly from its standard BPS baby counterpart. They are two rather different theories 
and their topological solutions do not have much in common. The main difference is that 
the vector BPS baby model in the massless limit does not possess finite energy solitonic 
solutions of the sourceless equations of motion. 
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B. Massive case 



Now we allow for a nonzero mass of the vector mesons. Then, only one equation of motion 
is modified. Namely, after the insertion of the ansatz into the static equations we find 



n\d T 



1 



1 + 



M 2 u - d x (2xuj' x 



(38) 



while the second equation obtained by variation with respect to the complex scalar remains 
unchanged 



/ 



f a/j, 



a-l 



0. 



(39) 



1. Compacton - a = 1 
In this case, the last equation can be completely integrated with the following result 

U=^r(X- X ) (40) 

nX 

for x < X and oj = for x > X, where at the moment X is an arbitrary constant - the size 
of the compacton. Now, we insert it into the first equation 



d x 



1 + / 2 



/i 2 M 2 
n 2 \ 2 



■x + 



n 2 \ 2 



(l + M 2 X) 



(41) 



which can be easily solved. After imposing the boundary conditions we get 

' -^x 2 + ^(l + M 2 X) x<X 

1 x > X 

If (i-t) x ^ x 



1 + f 2 {x) 



u[x) 



x>X 



where 



X = j^(-l + ^Jl + 2M 2 X ) 



(42) 



(43) 



(44) 



In the limit of small mass M — > we find that X — > X and we re-derive the former result. 
One can also check that 

X < X (45) 

and the size of the compacton tends to zero as M — > oo. So, the massive vector meson 
squeezes the compact baby skyrmion. 

Finally, one can again prove that our solution solves the field equations with the same source 
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term as in the massless case. This implies that, again, the topological charge is screened and 
that the solution is not a genuine topological soliton, although in the massive vector meson 
case there is no general physical principle which would enforce this result. And, indeed, we 
shall see in the next section that genuine finite energy solitons may exist for massive vector 
mesons. 



2. Exponentially localized solution - a = 2 case 



In order to find solutions for other values of the parameter a we differentiate ( 138]) and then 
insert ( )39l) . This gives a second order ODE for the profile function 

a.— X\ 7i ^-9 / no \ a— 1 

(46) 



1 



d 2 , 
x 1 1 + f 2 



x 



P 



M 2 



P 



l + P 



2X \1 + P 



This equation is quite complicated and, in general, no analytical solution can be found. 
However, there is one exception, namely the case with a = 2. Then, we find 

P \ M 2 f 2 „ , . n 2 \ 2 



0. 



x 



1 + P X l + f\ 



2X 1 + p 



X (a = 2) 



4/x 2 



p 2 

This equation can be brought to a linear form after the substitution g = -jj-p and x 



(47) 

2JQ 



M 



y 



d 2 y [(l + Py)g)]=g, (3 



V2 



X M 



(48) 



or 



g"(l + 0y) + 2/3g' -g = (49) 

with the boundary conditions g(0) = 1 and g(oo) = (here we already anticipate that the 
solutions are not compactons). Then the solution is 

1 K x (fvTT 



g(y) 



(50) 



M2/0) VT+Wy 

where K\ is the modified Bessel function of the second type. Hence, g(y) is a smooth function 
which goes to zero exponentially at spatial infinity. The same happens for the derivative of 

^ 1 K.dVT+fSy-) 



oj.„ 



(51) 



Now, we want to describe the energy of this configuration. In the field g(x) and variable x, 
the energy can be written as 

1 



E 



2tt / dx 
Jo 

POO 

2ix I dx 
Jo 



ji 2 g 2 — xtu 2 — -M 2 co 2 — n\cog a 



2 2 

V 9 



-n\cjg x 



d 

dx 



lxuu x 



(52) 
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where we used a partial integration and the (non-dynamical) equation 

- -4- (2xu x ) + M 2 u = -n\g x (53) 
ax 

for uj to eliminate the non-dynamical field u>. The resulting term — (l/2)n\ujg x (where u 
must be expressed by its solution) corresponds to the Coulomb self energy in the massless 
vector meson case, so we might call it the "Yukawa self energy" in the case of a massive 
vector meson. The corresponding Yukawa self-energy density is, in fact, positive semi- 
definite, because u(x) > and g x (x) < 0. Finally, the total derivative term gives zero for 
massive vector mesons because xuu x is zero both at x = and at x = oo. 

The energy is, in fact, of the BPS type, that is, it may be expressed as a total derivative 
with the help of the static field equations. Indeed, using 

2/i 2 

Wx = r9 (54) 

n\ 

(which follows from Eq. (|39|) ) the energy may be re-expressed like 

E = 2tt— ^ dx(<4 + uu xx ) = -- - u(0)u x (0). (55) 

Further, g(0) = 1 leads to u x (0) = — (2fi 2 /n\), whereas from Eq. (153]) . which in terms of 
the independent variable y may be written like 

M 

g y {l + Py) + (3g = --j=-uj (56) 



we get 

w(0) 



( 2 m (|) J < 57 > 

where we used that the y derivative of g(y) is 

/ 2 yi^ \ / 2v ^ \ / 2 yiT^ \ 

ft = -g 1 fej _ Ao Lxj + K ' izzzl. ,58) 

Putting everything together, we get the following energy expression 

g = ^f_g + *°(*) + *» W| ( 
" 2A 'i(I) 



where 





, = = ( 61 ) 

a/X^M nAM v ; 
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and we performed a small (3 (i.e., large n) expansion in the last line. We find that the BPS 
energies are not exactly linear in the topological charge n. Instead, they have sub-leading 
corrections in inverse powers of n, and the leading correction (proportional to (3, i.e., n^ 1 ) 
has a negative sign. This implies that the masses of small solitons of charge n = 1 are 
less than \jn times the masses of large solitons with charge n, and the large solitons are 
therefore energetically unstable w.r.t. decay into smaller ones, see Fig. [U It might still 

E/n 

5 - 

3 - 



2 4 6 8 10 

FIG. 1: Energy per topological charge n, as a function of n, for a = 2, and for the parameter 
values M = u = A = 1. The sublinear corrections reducing the energies of small solitons are clearly 
visible. 



be possible that there exist higher charge solitons without rotational symmetry with lesser 
energy, which would then be less unstable or even stable. All this is quite different from 
the case of the standard BPS baby skyrmions, where already the rotationally symmetric 
solitons saturate a bound exactly linear in the topological charge and are, therefore, stable. 
We remark that only for the case a = 2 we were able to express the on-shell energy as a 
total derivative, which implies that the spherically symmetric soliton solutions saturate a 
BPS bound. We were not able to find a similar result for other potentials. 



3. Exponentially localized solution - a € (1,2) case 



For other values of a G (1, 2) the power expansion at the vacuum value / = or g = 
in (1461) also leads to exponentially localized solutions. This follows from the fact that for 
g — > the field equation 



(J 2 

v 



g + Pyg 



a-l 



can be approximated by 



yg 



a-l 



g 



-a-l 



a-l 



(62) 
(63) 
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as g < g a 1 when < (a — 1) < 1. Then, the asymptotic solution, again, has an exponential 
tail 

2_ fW 

g a-1 A/ 8 

9 — r- (64) 

Obviously, the vector meson field as well as the energy density are exponentially localized, 
as well. 

An example of such a solution is given by a = 3/2. In this case we solve numerically the 
equations by shooting from the center. In this way we find the solutions presented in Fig. 
[2J where the initial values (i.e, values of the free integration constants) for this solution are 

g '(0) = -3.201, w(0) = 0.92043 (65) 




FIG. 2: Solutions g,u and their derivatives for a = 3/2, for the parameter values M = \i = n = 
A = 1. 
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4- Power-like localized solutions - a > 2 



It is easy to find that for a > 2 we get a power-like tail for the profile function 



/(*) 



x 



(66) 



which is exactly equal to the behavior of the profile function in the original BPS baby Skyrme 
model 10J. Then, one can obtain the asymptotic behaviour of the meson function 



x 



ix^-i 
xj 



(67) 



This behaviour is confirmed by the numerical analysis. For example, for a = 3 and for the 
parameter values M = /i = A = n= lwe find for the free integration constants for a soliton 
solution 

g '(0) = -3.94465, w(0) = 0.8327, (68) 

which gives the solutions presented in Fig. [3j 



1.0 



0.5 



0.0 



-0.5 
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FIG. 3: Solutions g, u and their derivatives for a = 3, for the parameter values M = fi = n = X = l. 
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C. Infinite mass limit 

We have demonstrated that there are no topologically nontrivial solutions to the vector 
version of the BPS baby Skyrme model with the old baby potential, unless some source 
terms are added. Thus, the vector BPS baby Skyrme model with the old baby Skyrme 
potential is qualitatively very different from the original BPS baby Skyrme model. 
From a phenomenological point of view this is not surprising. The BPS baby Skyrme model 
is based on the assumption of very (infinitely) heavy "baby" pions in the sense of a complete 
suppression of linear "pion" excitations. Physically, however, pions are the lightest mesons. 
Hence, in order to maintain the same physical hierarchy of masses, in this limit also other 
mesons should be infinitely heavy. 

On the other hand, in the vector BPS baby Skyrme model the vector mesons possess a well 
defined finite mass M while we still have infinitely heavy baby pions (no sigma term for the 
topological field). Thus, both models refer to a distinct mass hierarchy for the underlying 
particle d.o.f. Obviously, the limit represented by the vector model is not a physical one. 
Therefore, the fact that we do not reproduce results from the BPS baby Skyrme model in 
its vector version is a nice observation. Otherwise, the BPS model could be viewed (by its 
vector version) as a model with a completely unphysical assumption for the masses of the 
baby mesonic particles. 

In fact, there is a limit which brings the vector baby BPS model back to the BPS Skyrme 
model i.e, all baby mesonic particles have infinite masses. We simply assume that M — > oo 
i.e., we omit the kinetic part for the w-mesons in the Lagrangian. Of course, we arrive 
at a model in which the vector mesons are no longer dynamical d.o.f but rather a kind of 
Lagrange multipliers 



1 
2 

Now, the field equations for the w-mesons are 



£ = -y(0 3 ) + i M 2 ul + x '°°^- ( 69 ) 



X ';B, (70) 



" M 2 ' 

and can be used to eliminate these fields from the Lagrangian. Finally we re-derive the BPS 
baby Skyrme model 

(7i) 

III. INTEGRABILITY AND CONSERVATION LAWS 

The last issue which needs to be investigated is how the interaction with the vector mesons 
influences the integrability properties. In fact, we will show that the model is integrable in 
the sense of generalized integrability The first family of infinitely many conserved 

currents is given by the expression 

j? = iG(uu) (uTtf, - wx^) (72) 
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where G is an arbitrary function of the modulus of the complex scalar field and ir^ and 7r M 
denote the canonical momenta 

= § = tt - < "(T^¥ (73) 



dC 



iW** . (74) 



<9w^ q (1 + M 2 ) 2 ' 
The equations of motion can be written as 

9x = vv« - 2iAw^ , i y;„3 ( 75 ) 

(i + |wi 2 j d 

and 

0^ = _fv>u - 2iXuu a e^ - U ^ mr (76) 

(l + \u\ z Y 

Then, 

dtfG = iG ' (™V + ™v) - M7 v) + iG - ud^v) + iG (u% - u%) . (77) 

Using the identities which follow form the form of the momenta 

ttX = 0, 7fX = (78) 

and the field equations we get 

d^f G = i{G'uu + G) («% - u^) . (79) 

But 

(1 + M 2 ) 2 



= iXu a e a ^—^— = (80) 



and, therefore 

d,f G = 0. (81) 

Notice that we have used only the field equation for the complex scalar which is the same in 
the massless and massive w-meson models. Hence, the currents are conserved in both cases. 
For the massless model we can construct another infinite family of conserved currents. Now 
they have the form 

j* = H(uu)F IJH ,(uu v + uu v ) = HiuujF^iuu) (82) 

where H is an arbitrary function of the modulus of the complex scalar. We find 

= HF^d v {uu) + H'F^{uu)d v {uu) + H(& t F lu/ )&' '(««). (83) 

The first two terms vanish due to the contraction of the antisymmetric tensor F^ v with two 
symmetric ones. Using the field equation for the mesons we get for the third term 

d»fi = -iXHe™? UaU P (uu v + uu v ) = 0. (84) 
^ (1 + \u\ 2 ) 2 
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IV. SUMMARY AND CONCLUSIONS 



The vector BPS baby model is integrable (in the sense of generalized integrability 17] ) 
and to a far extent solvable for arbitrary values of the topological charge carried by a solitonic 
solution. So, in these aspects it reflects very well properties of the BPS baby Skyrme model. 
In the case of the old baby Skyrme potential, which may be viewed as a mass term for the 
Skyrme field, the BPS and vector models lead to quite different solutions. Namely, the vector 
model (massive as well as massless) requires a delta source term located at the boundary of 
the compacton, which equivalently can be stated as the appearance of a C-compacton, while 
in the original BPS baby model we get a C-compacton, that is, a solution to the sourceless 
equations of motion. So, strictly speaking, the vector BPS baby model with this potential 
does not support sourceless solitonic solution. This should be contrasted with the analogous 
situation for the full baby Skyrme model with the old baby Skyrme potential, where the 
soliton solutions of the original model and its vector meson version are very similar, see 



12]. As we pointed out already, the different behaviour in our case is a welcome result 
from the point of view of phenomenological applications to baryon physics. If the same 
happens in (3+1) dimensions (for the vector BPS Skyrme model) then the BPS Skyrme 
model would preserve the mass hierarchy of the mesonic fields in QCD. Otherwise, the BPS 
Skyrme model (infinitely heavy pion approximation) could be described by a model with 
finite massive vector mesons, and this is not acceptable from a phenomenological point of 
view as the pions are the lightest mesonic degrees of freedom. Undoubtedly, this should be 
checked in the future. 

For other potentials, where the "pion" potential cannot be interpreted as a mass term, we 
found that in some cases the soliton solutions in the original BPS baby Skyrme model and 
the massive vector model are more similar. Concretely, for a > 2 we observe exactly the 
same localization of solutions for the same values of the potential parameter a (exponential 
for a = 2, or power-like for a > 2), i.e., the same approach to the vacuum. For the specific 
case a = 2 we were able to find exact solutions and deduced from these solutions that, 
still, there remain some important differences. The solitons of the vector model for a = 2 
saturate a BPS bound, like the solutions of the BPS baby Skyrme model, but the nature 
of this BPS bound in the vector model is quite different. It is, in fact, a nonlinear function 
of the topological charge which constitutes a rather interesting result by itself (more recent 



results on nontrivial BPS bounds can be found, e.g. in [l8|, ll9|). As a consequence, a 
soliton with topological charge n is heavier than n solitons with charge one and is, therefore, 
energetically unstable against a decay into smaller solitons. This is completely different from 
the situation in the original BPS baby Skyrme model, where the solitons saturate a BPS 
bound which is linear in the topological charge and are, therefore, stable. The standard, 
inear BPS bound is, in fact, quite important for applications to hadronic or nuclear physics 



1, |3J, |20|, (21|. 



For a € (1, 2), the compactons of the original BPS baby model transform into exponentially 
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localized solitons with the tail becoming steeper and steeper as one approaches a = 1. It 
follows that the compactons, which exist in the original BPS baby Skyrme model for some 
potentials, completely disappear in the vector version. They either become solutions of the 
field equation with a source term (a = 1), or they receive exponential tails (a G (1, 2)). The 
disappearance of the compactons may perhaps be understood from the fact that the vector 
meson field possesses the standard kinetic and mass terms, which might prevent the solutions 
from becoming compactons. However, in the gauged BPS baby Skyrme model the addition 
of the standard kinetic term for the abelian gauge field does not lead to any problems with 



compactons [19j]. Therefore, the detailed form of the coupling of the additional (standard) 
field to the BPS baby model also plays a significant role. 

In the case of the massless model, all potentials with a > 1 generate large distance solutions 
with a logarithmically IR divergent energy, and this result can be understood easily by a 
comparison with the electrostatic Coulomb energy in two space dimensions. 
We summarize the types of topological solutions (baby skyrmions) we found in the three 
models in Fig. HJ 



BPS baby 



a=l a=2 



C'-compacton 

exponential power-like 



Massless vector 

a=l a=2 



C-compacton 



logarithmical 




C-compacton 



FIG. 4: Comparison of types of solutions in the BPS baby Skyrme model with its massless and 
massive vector counterparts. 
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